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Gamma-ray bursts: Potential sources of ultra high energy cosmic-rays 
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The arguments suggesting an association between tiie sources of cosmological gamma-ray bursts (GRBs) and 
the sources of ultra-high energy cosmic rays (UHECRs) are presented. Recent GRB and UHECR observations are 
shown to strengthen these arguments. Predictions of the GRB model for UHECR production, that may be tested 
with large area high energy cosmic-ray detectors which are either operating or under construction, are outlined. 



1. Introduction 

The origin of UHECR's is one of the most ex- 
citing open questions of high energy astrophysics 
|1I2| . The extreme energy of the highest energy 
events poses a challenge to models of particle 
acceleration. Since very few known astrophys- 
ical objects have characteristics indicating that 
they may allow acceleration of particles to the ob- 
served high energies the question of whether 
GRBs are possible UHECR sources is of great 
interest. Moreover, since the GRB model for 
UHECR production makes unique predictions, 
which differ from those of other models @] , the de- 
sign and analysis of future UHECR experiments, 
as well as of high energy neutrino and photon ex- 
periments, may be affected by the answer to this 
question. 

The phenomenology of GRBs, bursts of 0.1 
MeV-1 MeV photons lasting for a few seconds 
[S] , suggests that the observable effects are due to 
the dissipation of the kinetic energy of a cosmo- 
logically distant, relativistically expanding wind, 
a "fireball," whose primal cause is not yet known 
[HEIHl- Waxman Milgrom & Usov [HIl and Vi- 
etri have independently suggested that ultra- 
high energy, > 10"'^^ eV, cosmic rays may be pro- 
duced in GRB sources. The model suggested in 
ref. [HI was based on two arguments. First, it 
was shown that the constraints imposed on the 
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relativistic wind by the requirement that it pro- 
duces observed GRB characteristics are similar to 
the constraints imposed on such a wind by the re- 
quirement that it would allow proton acceleration 
to > 10^° eV. Second, the energy generation rate 
of 7-rays by GRBs was shown to be similar to the 
energy generation rate required to account for the 
observed UHECR flux ,9.12t . 

In §|2we briefly describe the model proposed in 
ref. for UHECR production in GRB fireballs. 
In §121 we discuss the implications of recent GRB 
and UHECR observations to this model (for a 
more detailed discussion, see JSl)- In §E]we dis- 
cuss predictions of the model that may be tested 
with high energy cosmic-ray experiments. For 
predictions of the model that may be tested with 
high energy neutrino detectors see j4ll4ll5j and 
references therein. Our conclusions are discussed 
in §□ 



2. Production of UHECRs in GRBs 

2.1. Proton acceleration 

General phcnomcnological considerations, 
based on 7-ray observations, indicate that, re- 
gardless of the nature of the underlying sources, 
GRB's are produced by the dissipation of the 
kinetic energy of a relativistic expanding fire- 
ball. A compact source, tq ^ 10^ cm, produces 
a wind, characterized by an average luminosity 
L ^ lO^^erg s^-^ and mass loss rate M. At small 
radius, the wind bulk Lorentz factor, F, grows lin- 
early with radius, until most of the wind energy 
is converted to kinetic energy and F saturates at 
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r ~ LjM(? ~ 300. Variability of the source on a 
time scale At ~ 10 ms, resulting in fluctuations 
in the wind bulk Lorentz factor F on a similar 
time scale, results in internal shocks in the ejecta 
at a radius r ~ ~ F^cAi ^ tq. It is assumed 
that internal shocks reconvert a substantial part 
of the kinetic energy to internal energy, which 
is then radiated as 7-rays by synchrotron and 
inverse- Compton radiation of shock-accelerated 
electrons. At a later stage, the shock wave driven 
into the surrounding medium by the expand- 
ing fireball ejecta leads to the emission of the 
lower-energy afterglow. 

The observed radiation is produced, both dur- 
ing the GRB and the afterglow, by synchrotron 
emission of shock accelerated electrons. In the re- 
gion where electrons are accelerated, protons are 
also expected to be shock accelerated. Thus, it is 
likely that protons, as well as electrons, are accel- 
erated to high energy within GRB fireballs. 

The internal shocks within the expanding wind 
are expected to be mildly relativistic in the wind 
rest frame, due to the fact that the allowed range 
of Lorentz factor fluctuations within the wind is 
from few x 10^ (the lower limit required to avoid 
large optical depth) to few xlO^ (the maximum 
Lorentz factor to which shell acceleration by ra- 
diation pressure is possible, e.g. |H]). This im- 
plies that the Lorentz factors associated with the 
relative velocities are not very large. Since inter- 
nal shocks are mildly relativistic, we expect our 
understanding of non-relativistic shock accelera- 
tion to apply to the acceleration of protons in 
these shocks. In particular, the predicted energy 
distribution of accelerated protons is expected to 
be dup/dtp oc e~^, similar to the electron energy 
spectrum inferred from the observed photon spec- 
trum. 

A power law energy spectrum with index 2 
has been observed for non-relativistic shocks (see, 
e.g., ref. ^Hl and references therein) and for rela- 
tivistic shocks It is consistent with that ex- 
pected theoretically for Fermi acceleration in col- 
lisionless shocks |1()I18| . although a first principles 
understanding of the process is not yet available 
(see, e.g. ref. ^Hl for a discussion of alternative 
shock acceleration processes). 

Several constraints must be satisfied by wind 



parameters in order to allow proton acceleration 
to high energy Cp. We summarize below these 
constraints. The reader is referred to refs. |9I4| 
for a detailed derivation. The requirement that 
the acceleration time, assumed to be compara- 
ble to the Larmor gyration time of the acceler- 
ated particle, be smaller than the wind expan- 
sion time, which also implies that the proton is 
confined to the acceleration region over the re- 
quired time, sets a lower limit to the strength of 
the wind magnetic field. This may be expressed 
as a lower limit to the ratio of magnetic field to 
electron energy density (2], 

ub/u, > OmTl,el,oL-',2. (1) 

where ep = lO^Oep^so eV, F = lO^-^Fj.s and = 
10^^L-y^52erg/s is the wind 7-ray luminosity. A 
second constraint is imposed by the requirement 
that the proton acceleration time be smaller than 
the proton energy loss time, which is dominated 
by synchrotron emission. This sets an upper limit 
to the magnetic field strength, which in turn sets 
a lower limit to F 

F > 13Qe'J^,%At-_l/\ (2) 

Here, At — 10^^At_2 s. As explained in the 
constraints Eq. and Eq. ((SJ hold regardless 
of whether the fireball is a sphere or a narrow 
jet (as long as the jet opening angle is > 1/F). 
The luminosity in Eq. ^ is the " isotropic equiv- 
alent luminosity" , i.e. the luminosity under the 
assumption of isotropic emission. 

Internal shocks within the wind take place at a 
radius « F^cAt. The constraint of Eq. is 
independent of At, i.e. independent of the inter- 
nal collision radius, while the constraint of Eq. Q 
sets a lower limit to the collision radius for a given 
At. This implies that protons may be accelerated 
to > 10^*^ eV regardless of the value of At, which 
may range from the dynamical time of the source 
(At ^ 1 ms) to the wind duration (At 1 s), pro- 
vided the magnetization and Lorentz factor are 
sufficiently large, following Eqs.n]and|5] 

At large radii the external medium affects fire- 
ball evolution, and a "reverse shock" is driven 
backward into the fireball ejecta and decelerates 
it. For typical GRB fireball parameters this shock 
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is also mildly relativistic (e.g. and its pa- 

rameters are similar to those of an internal shock 
with At ^ 10 s. Protons may therefore be accel- 
erated to > 10^*^ eV not only in the internal wind 
shocks, but also in the reverse shock 21 4 . This 
implies that proton acceleration to > lO^*' eV is 
possible, provided the constraints of Eqs.n]and|21 
(with At 10 s) are satisfied, also in (the cur- 
rently less favorable) scenario where GRB 7-rays 
are produced in the shock driven by the fireball 
into the surrounding gas, rather than by internal 
collisions, as suggested, e.g., in [221 ■ 

The constraints given by Eqs. and ( 
are remarkably similar to those inferred from 7- 
ray observations, based on independent physical 
arguments: F > 300 is implied by the 7-ray 
spectrum by the requirement to avoid high pair- 
production optical depth, and magnetic field close 
to equipartition, us/ug 0.1, is required in order 
to account for the observed 7-ray emission 6 7 8. 
This was the basis for the association of GRB's 
and UHECR's suggested in 0. 

It should be noted here, that the only assump- 
tion related to the acceleration process made in 
the derivation of the constraints ^ and Q is 
that the acceleration time is comparable to the 
Larmor gyration time of the accelerated particle. 
These constraints are valid therefore not only for 
Fermi-type shock acceleration, but rather to any 
electromagnetic acceleration process. This point 
is illustrated, for example, by the results of the 
analysis of ref. |23| , where a non Fermi- type elec- 
tromagnetic acceleration process is considered. 

2.2. gamma-ray and UHECR energy pro- 
duction rates 

The GRB model for UHECR production was 
suggested prior to the detection of GRB after- 
glows [55. Estimates of the rate of GRBs were 
based at that time on the 7-ray flux distribution, 
and ranged from ~ 3/Gpc'^yr to ~ 30/Gpc^yr 
|26| . The estimated average 7-ray energy release 
in a single GRB, based on a characteristic peak 
flux of ~ lO^^erg/s was ^ 10^^ erg. These es- 
timates were subject to large uncertainties, since 
the 7-ray luminosity function as well as the evo- 
lution of GRB rate with rcdshift were poorly con- 
strained. Based on the rate and energy estimates. 



the rate of 7-ray energy generation by GRBs was 
estimated to be ^ lO'^^erg/Mpc^yr. This rate is 
similar to the energy generation rate in cosmic- 
rays of energy in the range of 10^^ eV to 10^^ eV, 
4.5 ± 1.5 X lO'^^erg/Mpc^yr, inferred from Fly's 
Eye and AGASA measurements of the UHECR 
flux available at the time the GRB model for 
UHECR production was proposed |12| . 

The determination of GRB redshifts, which was 
made possible by the detection of afterglows, al- 
lows a more reliable estimate of the GRB energy 
production rate, and the increased exposure of 
UHECR experiments, provided by AGASA and 
HiRes, allows a more accurate estimate of the 
UHECR energy production rate. The implica- 
tions of these more recent observations is dis- 
cussed in § 13.21 

3. Implications of recent GRB and 
UHECR observations 

3.1. Proton acceleration 

Afterglow observations 24 lead to the confir- 
mation of the cosmological origin of GRBs and 
confirmed standard model predictions of after- 
glow that results from synchrotron emission of 
electrons accelerated to high energy in the highly 
relativistic shock driven by the fireball into its 
surrounding gas 6 7 8. Afterglow observations 
provide therefore strong support for the underly- 
ing fireball scenario. In addition, afterglow ob- 
servations provide important information on the 
values of model parameters that enter the con- 
straints given by Eqs. and 

Prior to the detection of afterglows, it was com- 
monly assumed that the farthest observed GRB's 
lie at redshift z ^ 1 |26I25| . Based on afterglow 
redshift determinations, we now know that de- 
tected GRB's typically lie at farther distances 
(e.g. "ST"). This implies that the character- 
istic GRB luminosity is higher by an order of 
magnitude compared to pre-afterglow estimates, 
« lO^^erg/s instead of w lO^^erg/s. This 
relaxes the constraint on magnetic field energy 
fraction given by Eq. 

In several cases, fast follow up afterglow obser- 
vations allowed the detection of radio and optical 
emission from the reverse shock (e.g. |28I29) '). 
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These observations provide direct information on 
the plasma conditions in the reverse shock, where 
acceleration of protons to high energy may take 
place. Two major conclusions were drawn from 
the analysis of the early optical and radio reverse 
shock emission. First, lower limits to the initial 
fireball Lorentz factors were inferred, in the range 
of r > 100 to r > 1000 . Second, the mag- 

netic field in the reverse shock was inferred to be 
close to equipartition, that is us/ue was inferred 
to be of order unity 30 28 . Early afterglow ob- 
servations provide therefore constraints on T and 
on us/ue which are (i) Independent of the con- 
straints derived from 7-ray observations; (ii) Con- 
sistent with the 7-ray constraints; and (iii) Are 
remarkably similar to the constraints of Eqs. 
and (0) , that need to be satisfied in order to allow 
proton acceleration to > 10^" eV. 

3.2. gamma-ray and UHECR energy pro- 
duction rates 
3.2.1. gamma-ray production rate 

Most of the GRBs are observed from z > 1, 
since they can be detected out to large redshift. 
This implies that the GRB rate density at z > 1 
is better constrained by the observations than 
the local, z — rate. The inferred local rate 
depends on the assumed redshift evolution. It 
is now commonly believed that the GRB rate 
evolves with redshift following the star-formation 
rate, based on the association of GRBs with type 
Ib/c supernovae. This association is based on 
the temporal and angular coincidence of several 
GRBs and type Ib/c supernovae |Hll32l33j . and 
on evidence for optical supernovae emission in 
several GRB afterglows Adopting the as- 

sumption, that the GRB rate follows the redshift 
evolution of the star formation rate, the local 
{z ^ 0) GRB rate density is -Rgrb(^ = 0) « 
0.5 X lO^^Mpc"^ yr~^. Given the current (sys- 
tematic uncertainties in the redshift) data, this 
rate is accurate to within a factor of a few j33| . 

The local energy generation rate in 7-rays by 
GRB's, e-y, is given by the product of Rgkb[z — 
0) and the average 7-ray energy release in a 
single GRB, e-y. [211 provide for 27 bursts 
with known redshifts, in a standard rest-frame 
bandpass, 0.02 MeV to 2 MeV. The average is 



= 2.9 X 10^'^ erg, with estimated uncertainty, 
due to the correction to a fixed rest-frame band- 
pass, of ~ 20% for individual bursts (and much 
smaller for the average). In calculating from 
this value of e^, the following point should be 
taken into account. e~f is the average energy for 
bursts with known redshift, most of which were 
localized by the BeppoSAX satellite. Since Bep- 
poSAX has a higher detection flux threshold than 
BATSE it is sensitive to « 70% of the 

bursts detectable by BATSE, for which the GRB 
rate Rg^b{z — 0) was inferred. Thus, the en- 
ergy generation rate by bursts detectable by Bep- 
poSAX is 

^7[0.02McV,2McV] ~ 0.7i?GRB(z = 0)e^ 

= lO^'^erg Mpc~3 yr~^ (3) 

The energy observed in 7-rays in the range of 
BATSE is dominated by photons in the energy 
range of 0.1 MeV to 1 MeV, and therefore reflects 
the energy of accelerated electrons over roughly 
half a decade of electron energies (recall that for 
synchrotron and inverse-Compton emission the 
photon energy is proportional to the square of the 
electron energy). Thus, the rate of energy gener- 
ation per logarithmic interval of electron energy 
is 

2 dTl^ . GRB 

~ ^7[0.02MoV,2McV] 

« lO-^^erg Mpc"^ yr-\ (4) 

The main uncertainty in determining e-y is re- 
lated to the uncertainty in the local GRB rate, 
due to which the value given in Eq.|21is accurate 
to within a factor of few. It should be noted that 
the numbers quoted above for the GRB rate den- 
sity and 7-ray energy release are based on the as- 
sumption that GRB 7-ray emission is isotropic. 
If, as now commonly believed, the emission is 
confined to a solid angle Af2 < 47r, then the 
GRB rate is increased by a factor (Ail/47r)^^ 
and the GRB energy is decreased by the same 
factor. However, their product, the energy gen- 
eration rate, is independent of the solid angle of 
emission. 
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3.2.2. UHECR production rate 

The cosmic-ray spectrum flattens at ^ 10^^ eV 
|37l38j . There are indications that the spectral 
change is correlated with a change in composi- 
tion, from heavy to hght nuclei j37.,39 40 (see, 
however, These characteristics suggest that 

the cosmic ray flux is dominated at energies < 
10^^ eV by a Galactic component of heavy nu- 
clei, and at UHE by an extra-Galactic source 
of protons. Also, both the AGASA and Fly's 
Eye experiments report an enhancement of the 
cosmic-ray flux near the Galactic disk at energies 
< 10^^-^ eV, but not at higher energies j42l43j . 
Fly's Eye stereo spectrum is well fitted in the 
energy range 10^^'^ eV to 10^^'^ eV by a sum 
of two power laws: A steeper component, with 
differential number spectrum J cx £~3-50^ domi- 
nating at lower energy, and a shallower compo- 
nent, J cx e~'^'^^ , dominating at higher energy, 
e > 10^^ eV. The data are consistent with the 
steeper component being composed of heavy nu- 
clei primaries, and the lighter one being composed 
of proton primaries. 

The observed UHECR flux and spectrum may 
be accounted for by a two component model |44) . 
with a Galactic component given by the Fly's Eye 
flt,^ oc £-3.50^ g^j-^j extra-Galactic compo- 
nent characterized by a local {z = 0) energy gen- 
eration rate of extra-Galactic protons of 

= 0.65 X lO^'^erg Mpc^^ yr-i. (5) 

Uncertainties in the absolute energy calibration 
of the experiments lead to uncertainty of w 20% 
in this rate 1 I . The spectral index, 2, is that 
expected for acceleration in sub-relativistic colli- 
sionless shocks in general, and in particular for 
the GRB model discussed in § 12.11 The model 
used in 0^ is similar to that proposed in jl2| . 
The improved constraints on UHECR spectrum 
and flux provided by the recent observations of 
HiRes do not change the estimates given in [T2] 
for the energy generation rate and spectrum, 
Eq. jSJ, but reduce the uncertainties. 

Comparing eqs. ^ and ^ we find that the 
observed UHECR flux and spectrum may be ac- 
counted for if GRBs produce high energy elec- 
trons and protons with similar spectra and rates. 



It should be emphasized here, that the ratio be- 
tween the energy carried by relativistic electrons 
and protons in collisionless shocks is not known 
from basic principles. Moreover, as explained 
above, the estimate of the GRB gamma-ray en- 
ergy generation rate is uncertain by a factor of 
a few. Thus, an exact match between the de- 
rived 7-ray and UHECR generation rates should 
not necessarily be expected. Given current un- 
certainties, the two rates should only be expected 
to be of the same order of magnitude. 

Finally, the following point should be men- 
tioned. We have so far discussed the contribu- 
tion of extra-Galactic GRBs to the UHECR flux. 
GRBs exploding within our Galaxy may also con- 
tribute to the Galactic cosmic-ray flux at lower 
energies, > 10^^ eV jlHl- However, in order for 
GRBs to make a signiflcant contribution to the 
flux at 10^^ eV, they should produce more en- 
ergy in ^ 10"'^^ eV protons than they are inferred 
(based on their gamma-ray emission) to produce 
in high energy electrons |46|. 

4. Predictions for high energy cosmic-ray 
experiments 

The initial proton energy, necessary to have an 
observed energy ep, increases with source distance 
due to propagation energy losses. The rapid in- 
crease of the initial energy after it exceeds, due 
to electron-positron production, the threshold for 
pion production effectively introduces a cutoff 
distance, Dc(ep), beyond which sources do not 
contribute to the flux above ep. The function 
Dci^p) is shown in Fig. (adapted from 02]). 
Since Dc{ep) is a decreasing function of Cp, for a 
given number density of sources there is a critical 
energy Cc, above which only one source (on aver- 
age) contributes to the flux. In the GRB model ec 
depends on the product of the burst rate Rgrb 
and the time delay. The number of sources con- 
tributing, on average, to the flux at energy tp is 

m 

Ajr 

N{ep) = —RGRBD,{ep)\[ep,D,{ep)] , (6) 

where T(ep, D), the spread in arrival time of pro- 
tons of energy ep produced by a source at dis- 



6 




Figure 1. Results of a Monte-Carlo realization 
of the bursting sources model, with ec — 1.4 x 
10^° eV: Thick solid line- overall spectrum in 
the realization; Thin solid line- average spectrum, 
this curve also gives Dc(ep); Dotted lines- spectra 
of brightest sources at different energies. 



tance D due to deflection by inter-galactic mag- 
netic fields, is pTz] 

T{ep,D)^ 10' e;',,DloaXMpcB% yr. (7) 

Here D = lOOZ^iooMpc, A = IAmpc Mpc is the 
characteristic length scale over which the field di- 
rection varies, and B — lQ~^Bs G is the char- 
acteristic field amplitude. The average intensity 
resulting from all sources is 

J{ep) = ^RGRBp^D.iep) , (8) 
47r d€p 

where dup/dsp is the number per unit energy of 
protons produced on average by a single burst 
(this is the formal definition of Dc{ep)). 

The critical energy Cc, beyond which a single 
source contributes on average to the flux, 
Air 

-^RGR.BDc{ecfT [ec, D,{ec)] = 1 , 
5 



(9) 



depends on the unknown properties of the inter- 
galactic magnetic field, r oc B'^X. However, the 
rapid decrease of -Dc(ep) with energy near lO^^eV 
implies that ec is only weakly dependent on the 
value of B'^X. In The GRB model, the prod- 
uct Rgrbt{D = 100Mpc,ep = lO^OeV) is ap- 
proximately limited to the range 10~^ Mpc~^ 
to 10^2 Mpc"^ (The lower hmit is set by the 
requirement that at least a few GRB sources 
be present at I? < 100 Mpc, and the upper 
limit by the Faraday rotation bound BX^/^ < 
10-**GMpci/2 g^^ g^^^ ^^^^ < 1/ GpcVr)- 

The corresponding range of values of is 
lO^^eV < Cc < 4 X 102"eV. 

Fig. n presents the flux obtained in one re- 
alization of a Monte-Carlo simulation described 
by Miralda-Escude & Waxman 47 of the total 
number of UHECRs received from GRBs at some 
fixed time. For each realization the positions (dis- 
tances from Earth) and times at which cosmo- 
logical GRBs occurred were randomly drawn, as- 
suming an intrinsic proton generation spectrum 
dup/dep oc e~^, and Ec = 1.4 x 10^°eV. Most of 
the realizations gave an overall spectrum similar 
to that obtained in the realization of Fig. ^when 
the brightest source of this realization (dominat- 
ing at 10^°eV) is not included. At ep < ec, the 
number of sources contributing to the flux is very 
large, and the overall UHECR flux received at 
any given time is near the average (the average 
flux is that obtained when the UHECR emissivity 
is spatially uniform and time independent). At 
ep > ec, the flux will generally be much lower than 
the average, because there will be no burst within 
a distance Dc{ep) having taken place sufficiently 
recently. There is, however, a significant prob- 
ability to observe one source with a fiux higher 
than the average. A source similar to the bright- 
est one in Fig. ^ appears ~ 5% of the time. 

At any fixed time a given burst is observed in 
UHECRs only over a narrow range of energy, be- 
cause if a burst is currently observed at some 
energy ep then UHECRs of much lower energy 
from this burst have not yet arrived, while higher 
energy UHECRs reached us mostly in the past. 
As mentioned above, for energies above the pion 
production threshold, ep ~ 5 x lO^^eV, the dis- 
persion in arrival times of UHECRs with fixed 
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observed energy is comparable to the average de- 
lay at that energy. This implies that the spec- 
tral width Aep of the source at a given time is 
of order the average observed energy, Ae^ ep. 
Thus, bursting UHECR sources should have nar- 
rowly peaked energy spectra, and the brightest 
sources should be different at different energies. 
For steady state sources, on the other hand, the 
brightest source at high energies should also be 
the brightest one at low energies, its fractional 
contribution to the overall flux decreasing to low 
energy only as Ddep)^^- A detailed numerical 
analysis of the time dependent energy spectrum 
of bursting sources is given in ,49.511) . 

5. Conclusions 

The main constraints that a relativistic wind 
(fireball) need to satisfy to allow proton accel- 
eration to > 10^° eV are given by Eqs. ^ and 
(|2Jl : The magnetic field energy density u b should 
exceed a few percent of the relativistic electron 
energy density Ue, and the wind Lorentz factor 
r should exceed « f 0^. As explained in § 12.11 
these constraints are independent of the details 
of the acceleration process, and are valid for any 
electromagnetic acceleration process. The sim- 
ilarity of these constraints and the constraints 
imposed on wind parameters, based on indepen- 
dent physical considerations, by 7-ray observa- 
tions were the basis for the association of GRB 
and UHECR sources suggested in 0. We have 
pointed out in § 13.11 that afterglow observations 
of GRBs imply a higher characteristic GRB lu- 
minosity than estimated based on 7-ray observa- 
tions alone, thus relaxing the constraint of Eq. Q 
on UB/ug. Moreover, early optical and radio af- 
terglow observations provide new constraints on 
wind parameters, implying large Lorentz factors, 
r > 10^ to r > 10^, and large magnetic field 
energy density in the fireball plasma, us/ue 1- 
These constraints are consistent with those previ- 
ously inferred from 7-ray observations, and with 
the constraints imposed by the requirement to al- 
low proton acceleration to > 10^° eV. 

We have also pointed out, in S 12.21 that GRBs 
produce high energy electrons at a rate (and with 
a spectrum) similar to the rate (and spectrum) at 



which high energy protons should be produced in 
order to account for the UHECR spectrum and 
flux. Thus, the observed UHECR flux and spec- 
trum may be accounted for if GRBs produce high 
energy electrons and protons at a similar rate and 
spectrum. 

In §21 we have pointed out that the local GRB 
rate implies that the rate of GRBs out to a dis- 
tance from which most protons of energy exceed- 
ing 10^° eV originate, ~ 90 Mpc (see fig. is 
^ 10~'^/yr. The number of GRBs contributing 
to the observed flux at any given time, eq. © is 
given by the product of this rate and the spread 
in arrival time of protons, due to the combined 
effect of stochastic propagation energy loss and 
deflection by magnetic fields. This time spread, 
eq. 0, may be as large as 10'' yr for lO'^" eV 
originating at 90 Mpc distance, implying that the 
number of GRBs contributing to the > 10^" eV 
flux at any given time may reach 10*. The 
upper limit on the strength of the inter-galactic 
magnetic field, combined with the low local rate 
of GRB's, leads to unique predictions of the GRB 
model for UHECR production |47I51| . that may 
be tested with operating |52) . under-construction 
and planned large area UHECR detec- 
tors. In particular, a critical energy is predicted 
to exist, lO^OeV < < 4, x lO^^eV, above 
which a few sources produce most of the UHECR 
flux, and the observed spectra of these sources is 
predicted to be narrow, Ae/e ~ 1: The bright 
sources at high energy should be absent in UHE- 
CRs of much lower energy, since particles take 
longer to arrive the lower their energy. 
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